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ABSTRACT 


The constant lift rotor (CLR) enploys a control input of 
pitch moment to several airfoil sections which are free to pivot 
on a continuous spar, allowing them to change their pitch to ob- 
tain the desired lift. A flap-lag-torsion flutter analysis of 
a constant lift rotor blade in hover was developed. The blade 
model assumes rigid body flap and lead-lag motions at the root 
hinge and each strip undergoes an independent torsional motion. 
The results are presented in terms of root locus plots of complex 
eigenvalues as a function of thrust. The effects of severed par- 
ameters (including structural danping, center of gravity and elas- 
tic axis offset from aerodynamic center, compressibility, pitch- 
lag and pitch- flap coupling) on the blade dynamics are examined. 
With a suitable combination of lag damper and pitch-flap coupling, 
it is possible to design a constant lift rotor blade free frern 
flutter instability. 
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Blade section angle of attack 
Angular deflections 
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Blade lock number, p£ r" 
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Rotor inflow ratio 

Nonrotating natural frequencies of blade 

Nondimensional coordinate, r^/i 

Mr density 
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H7TRDDUCTI0N 


In conventional rotors, lift variation is achieved by pitch 
control input to blades. For a constant lift rotor (CLR) , the 
lift variation is achieved through a control input of pitching 
moment to the blade and allowing the blade to change its pitch. 

This helps to alleviate the oscillatory loads on the rotor. 

Constant lift rotor employs a finite number of segments pivot- 
ally mounted on a continuous spar and the pitch of each segment is 
determined by the balance of centrifugal, aerodynamic, control, and 
frictional forces. Each one of the airfoil strips is directed 
through a control rod to achieve a desired amount of lift. There 
is chordwise offset of elastic axis from aerodynr.ic center (forward 
direction) and the strips can float freely torsionally. 

In the present paper, the flap-lag-torsion flutter of a con- 
stant lift rotox blade in hovering is investigated. The equations 
of motion for the shaft-fixed dynamics are derived for a blade with 
finite number of spanwise strips. These equations are linearized 
about a trim static solution in hover. The blade is assumed to have 
rigid body flap and lead- lag motions at the root hinge and, also, 
each one of the segments can undergo independent torsional motion. 
Qu$i-static airfoil characteristics are used to obtain aerodynamic 
forces (stalling is not considered, however) . The results are pre- 
sented in form of root locus plots of complex eigenvalues as a func- 
tion of c^/or The effects of structural daitping, eg and elastic 
axis offset from aerodynamic center, compressibility correction, 
pitch- flap and pitch-lag coupling are studied. 
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EQUATIONS OF MOTION 


Fig. 1 presents the blade configuration considered for the 
analysis. The blade consists of N rigid strips, M connected 
to spar through torsional springs, and its flap and lead- lag 
stiffnesses are represented by springs at hinges offset by a dis- 
tance el from the hub. The hinge sequence is flap inboard, lead- 
lag and then torsional motions outboard. The flap angle $ is pos- 
itive up, the lead-lag c is positive forward, and the pitch angles 
0j_ are positive nose up . The hub has a precone angle 8p. The 
blade rotates at a constant rotational speed ft. 

The linear equations of motion are derived by assuming the 
response consists of small perturbation motion (8,S,6]_, . . .0 N ) 
about a steady deflection (8 0 /So) . In general, terms up to second 
order are retained in flap and lag equations, terms up to third or- 
der are retained in torsion equations. 

Flap eqn. : 


N 


6 + 2(6 0 +6 p )C + (1+ f etv 8 )8 + g e 6 - Z R^W 

-If - w 2 


(1) 


Lag eqn. : 


N 


C - 2(8 0+ 8 p >8 + (|e+v ? 2 )5 +■ - I 'V.<6 o +6 p )0 i 


N 


+ 4 z X s. (£.+£. ,,){0 .0.+2(B +8 )0-} = Q /U 
2 £ *cg^ i' s i ^l+l' oi i o p i CD 


( 2 ) 


Torsion eqn. : (i^ strip) 


VW +(1+ V )e i* < 6 o +e p )! ' c o s+ 9e. 5} 


+ ? * 03 . ( *cg. Vi" I V«i + W <6 +6,+ I s i< f M +5 i> [6 oi ! ’- 2 
1 1 


= W 


(3) 
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where is the ratio of torsional inertia of i^ strip to blade 
inertia 1^; is the ratio of the chordwise offset of og and 
pitch axis to rotor radius (positive toward the trailing edge) ; 


g., are the structural dartping coefficients; v 0 , v_, v 0 

P ? P 5 

are the nonrotating natural frequencies of the blade (divided by fl) . 
It is assured in the above derivations that the inertia properties 
are uniform within each strip and ^ and s^ respectively repre- 
sent strip beginning from hinge and strip width in terms of blade 

length i. The Q_, Q and Q. are the perturbation aerodynamic 
p K 

moments. 

The trim equations are: 


3_ . 2 


2 3 


N 


(1+ |efv 6 Oe o + (1+ |e)6 p - Vg e s - (8 0 +8 p ) J | x^ . Si (£ ± +{ i+1 ) e oi 


= % /V 


,3.. 2,. 2_ 3 


N 


<2^ >«o - 5 s " 1 1 <e o + V 'o 


(4) 


where 6 and c are the initial settings of flap and lag hinge 
s s 

springs (relative to the preconed hub) . Q 0 and Q_ are the 

o C o 

steady aerodynamic moments at the hinges. 

The aerodynamic forces are obtained using quasi-steady air- 
foil theory. The section lift, drag and moment about aerodynamic 
center are: 
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JsPT 2 cc £ (a,M) 

*sPV 2 cc d («,M) 

hPV*c 2 c 

m 


( 5 ) 


Hie section lift and drag coefficients c £ and c^ are functions 
of section angle of attack and Mach number. 


C £ = (c Q +c 1 ot)CR 

c d = 6 0 + 6 i a + 6 2 a2 + Ac d 


Hie compressibility correction CR and Ac^ are 


where 



M < M 
c 



1-M 

1-M 

c 


M > M 
C 



(7) 


= 0 
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div 


div 


. 26 (M ,-M) 
a 


= 0 


M > M . 
a 


M c and are respectively lift divergence and drag divergence 

Mach numbers. 
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Fig. 2 shows the section aerodynamic environment. The flew 


velocity oenponents along the shaft axes are and and 
0 Q ^ in the pitch for the i^ strip. The resultant flew velocity 


is then V = ^ u 2 4U 2 . The resolved aerodynamic mcments in 

T + p 


shaft axes are 


N = L oos<|> - D sin$ 


C = L sin<|> + d oos<t> 


( 8 ) 


M = M 


ac 


where <p is the induced angle of attack, tan /il, 

r A 

The perturbation section aerodynamic forces and pitch moment 
are 


an - Jscc[6u^{- ( Vt -u^ d ) + (c^ ) - c d v 

oc oc M 

2 

-“(■y-V 

2 

+ an { ^ ^ -ujp d > + U -v~ ‘ c 1 #fc n M > + c ,y 

n a M 

Vt 

- (c d +Mo dH )) + «e 0 t v «*^r UpC^}}] 


U, 




AC = Hpc [«lpf- 7T- (UpC^ -»u^ d ) + (c^+Mc, ) +c 

oc a M 

+ t (Vi + '¥’d 1 + V 1 + c d v 

oc oc M 

. 2 


„v 


(9) 


+ 60 o {V(u p c £ + Vd )> 


(10) 
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+ 4 Wm + 2 Vm } 

* CC 

+ ^{V 2 ^ } (11) 


and the steady forces are 

N q = J 5 Pc{c ? U T V-C d U p V} 

C o = ^iV^dV* 


(12) 


The perturb 7 cion aerodynamic rrrxnents required for equations of 
notion are 


°6 “ A d-^ 0 2 )rdN - u o cN o - re .C o 

+ X A C c (dN oose ^-HK -taO ) + X A C(N cose ^ sine ^Mr 
e e 

Q - - ; £ {rdc f (13) 

s o 

Q e< = - i <x A oose oi dN+x A sine^dC+dM+M^Jdr 
i e e 


and the steady aerodynamic moments are 


Qg = ? 0 (N 0 cxxse^ sine oi )}dr 

00 e 

Q_ = / £ rCdr 
o 


(14) 


where r is the radial distance from hinge and x h 


is the chord- 


wise offset of pitch axis from aerodynamic center (positive for- 
ward) . is the noncirulatory aerodynamic pitch mcxnent due 

to unsteady thin airfoil theory. 


■fa =in t , S 1 2 o 3 (r(J + -^)S-r(| + -^)e 

X A X» 2 

- c (- 3 l + if + 


( 14 ) 


The flow velocity oonponents for hover are 


Steady: * T = «£{ (C+e) - ^(^Bp) > 


U p = lull + £(%w p )e 0 > 

Perturbations: 6^ = S2H{CC - 5(B 0 +B p )B} 


&Up = OAUB + €(8 0 +3 p )C + 5C 0 6> 


and perturbation pitch for small angle assumption 


(15) 


(16) 


( j sch-x a ) 

60 = 0 . + — 0 . 

° 1 (C+e)Jt 1 


(17) 


where £ = r/l and A^ is the wake induced inflow ratio (divided 
by fli.) and is assumed to be uniform within each strip for hover- 
ing condition. The A^ is obtained by equating the thrust for 
i strip obtained from momentum theory and blade element theory 
for a specified load distribution: 
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( 18 ) 



where c^, is the thrust coefficient and f^ is the fractional 
thrust produced by 1 th strip (AT\ = f \T ) . The pitch settings 
are obtained as 



The resulting equations of motion for N-strips blade are 



The coefficients of these matrices sure defined in appendix. 
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FITCH-LAG AND PITCH-FLAP COUPLING 


Par constant lift rotors, control input to individual blade 
strip is the pitch moment; thus, it is assumed here that the pitch- 
lag and pitch-flap coupling moments for each one of the strips are 
proportional to their torsional stiffness. 

FOr 1 1 strip 



" R 3 K i+2 i+2 


3 + 


R ? K i+2 i+2 C 


The pitch-flap coupling R is positive flap up, resulting in 
pitch down control mcment P and the pitch-lag R is positive lag 
back, resulting in pitch down control moment. ^The inclusion of 
these coupling moments modifies the stiffness matrix of resulting 
equations. 

FOr i^ 1 strip 

Pitch-flap coupling 


K i+2 1 (K i+2 + R B K i+2 i+2 

Pitch-lag coupling 

K i+2 2 = (K i+2 2 } 0 R ? K i+2 i+2 

where ( ) Q represents the blade without coupling. 
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RESULTS AND DISCUSSION 


The flutter stability is examined for a constant lift rotor 
blade with lock number y = 8.65? solidity ratio a = .088? pre- 
oone Bn = 3. 5° ? hinge offset e = .04? chord to radius ratio 
c/R = .04 and no droop or sweep. The following airfoil charac- 
teristics are used 


c i - 5 - 7a 

c. = .008 + .023 01 + .076 a 2 


The constant lift rotor blade consists of five strips with the 
following properties 


stiify 

1 

Z 

3 

H 

5 

stiifj width in 
fULdius , 

•6 

•/ 

■1 

• i 

•/ 

Tonxon intsLa./ FLu{> 
(a. Rmi 

•00018 

■00003 

■ 00003 

■00003 

.00003 

Thrust teJxOjf; 

■3S 

• IS 

•n 

■H 

•/*/ 


R>r ocnpressibility cojrrections , the lift divergence and drag diver- 
gence Mach numbers are assumed to be 0.7. Results are also presented 
for a conventional rotor blade (single strip) for ccmparision. Figure 
3 shows the trim solution, the collective pitch as a function of thrust 
for constant lift rotor as well as conventional rotor. For dynamic 
results, the nondimensional eigenvalues (real and imaginary) are plot- 
ted in the complex plane for increasing thrust c_/a. 


- 10 - 
















Figures 4-6 present the dynamic stability of a single strip 
blade. The first case considered is a <x>nventional blade with 
torsion frequency of wq = s/rev and with no elastic axis/aero- 
dynamic center offset = 0) (Fig. 4). The blade is stable 

except near zero thrust level where lead-lag mode gets into a very 
weak instability. However, with zero torsion frequency ( 0)9 - 0) 
and zero elastic axis offset (x* = 0 ) , the blade becomes very 
unstable (Fig. 5) . Hie lead-lag^instability expands up to higher 
thrust levels; and at still higher thrusts, torsion divergence takes 
place. With inclusion of elastic axis effect of XAe = -12c, the 
blade beccmes torsionally stiff (due to aerodynamic forces) (Fig. 

6 ) . The lead-lag mode is still unstable and it beccmes more and 
more violent with increasing level of thrist. The flap mode also 
gets into a weak instability at lower thrusts. 

Figure 7 shows the results for a constant lift rotor blade 
with five strips, freely floating torsionally (uq^ = 0 ) , and 
with elastic axis offset (x^g. = *12c). All the seven eigen- 
values are plotted here. The nature of too lowest daitped modes 
which happen to be flap and lag modes is very similar to single 
strip case (see Fig. 4). The torsion modes of different segments 
are quite stable and therefore in subsequent figures only flap and 
lag modes are plotted. 


the effect of elastic axis offset from 
J on the blade dynamics of constant 
5r Xp^ - .15c, the flap and lead-lag 

modes are hardly different from 1 those for the blade with 
- . 12 c. The offset Xp^ primarily affects the different 


In Figures 8 (a) - (b) , 
the aerodynamic center (x^ 
lift rotor is presented. I 




torsion modes, in fact, with increasing x^ the torsional modes 
become stiffer, as expected. 


Figures 9 (a) -(b) show the influence of the eg offest from 
elastic axis, Xcg (positive toward trailing edge) on the blade 
dynamics. For Xcg-^ = • 12c (eg coincidental with aerodynamic 

center) , the flap and lag inodes are changed very much and also 
torsional frequencies are reduced from that of the blade with 
X = 0. (Fig. 9(a)). Also, for sane lower thrusts, the lag 
°*i 

mode gets into static divergence. For x~_ = *06c (c.g. lying 

midway between the elastic axis and the aerodynamic center) , the 
nature of flap and lag modes is quite similar to that of blade 
with x^. =0. (Fig. 9(b)). 
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In Figures 10-11, the influence of structural danping on 
the dynamics of constant lift rotor blade is shown. The flap 
mode instability, which was mild, can be easily stabilized with 
a low level of flap damper (Fig. 10) . On the other hand, the 
lag mode instability becomes increasingly violent with higher 
thrusts, needs a fairly big lag damper to stabilize it (Figs. 

11(a) and 11(b)). 

Figures 12-18 shew the effects of pitch-flap and pitch-lag 
coupling on the blade stability. The inclusion of positive pitch- 
flap coupling (pitch down moment with flap up) increases the cross 
coupling stiffness (k. +2 ,) for I th strip. This destabilizes 
the flap mode somewhat 1 (Fig. 12 (a) ) . The opposite effect is 
seen with negative pitch-flap coupling which stabilizes the flap 
mode. (Fig. 12 (b) ) . The addition of positive pitch-lag coupling 
(pitch down moment with lag back) stiffens the lag mode slightly, 
however, with little effect on its stability. (Fig. 13 (a) ) . Again, 
with the negative pitch-lag coupling, the lag mode gets softened 
but the instability region is nearly the same (Fig. 13 (b) ) . 

Figure 14 shows the effect of oonpressibility for a blade 
with tip Mach number of .6 (x A = .12C). The general behavior 

e 

is similar to the blade reglecting the oonpressibility effects, 
particularly for low thrusts. 

Figure 15 presents the blade dynamics with a suitable combi- 
nation of lag damper (g^ = .5) and pitch-flap coupling (R^ = -.3) 

for a constant lift rotor (x. - .12C, x = 0). The blade is 

A eg. 

e i 1 

quite stable in the covered range of thrust. The same combination 
of lag damper and pitch-flap coupling is also used for another 
blade configuration (x. = .12C, x~. = .06C). In Figure 16 

A Cm • 

e i 

the blade is shown to be stable except at very high thrusts 

Op 

(— > .2) the lag mode gets into static divergence. 
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CONCH JSICNS 


The flap-lag-torsion flutter of a constant lift rotor in 
hover has been investigated. The CLR blade consists of a finite 
nurber of strips pivotally mounted on the spar and their torsional 
stiffness is attained through the elastic axis offset from the 
aerodynamic center. The perturbation equations of motion were 
derived, retaining the higher order steady terms. The dynamic 
results for multi-strip constant lift rotor are quire similar 
to those of a single strip blade under the same environment. 

The effects of severed, parameters on the blade dynamics 
were examined, including structural dairying, eg and elastic 
axis offset from aerodynamic center, compressibility correction, 
pitch-lag and pitch- flap coupling. With a suitable combination 
of lag danper and pitch-flap coupling, it is possible to design 
a constant lift rotor blade free from aeroelastic instability. 
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APPENDIX 


Inertia Matrix "M" 



= to 
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2. Damping Matrix "C" 
(a) Structure 


c " " ^ 
c , 2 = 2 f£ + £) 

c «/ f = 0 

C S) = - z d + ^) 
C 22 = 
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= 3 x ‘i. s ;(?i + ?u,)0 l+f) 
C fc, = 0 

c hi = -3* % s l (z + f.j(/ 2 +p p ) 

^ hh ~~ ^n t - %• 


(b) Aerodynamics 


+ %X <l R z C 0 S 6 oi }] 

c ,z =-!|f r, { e T % *> 3 * + r,{^ 4*^ 

+ X ^ COS6 oi}] 

^ih — f [ ^{tsC-i & + I 2 { eg s+2 i R 6 +x llt i e R s a)s(i']jl 
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3. Stiffness Matrix "K" 

(a) Structure 

k„ = m ^ 

k u = 0 

K >k « ~{ l 

k 3l — 0 

K iz = ie + if 
*zk = + 

K U = - \ i-f X % S c ( ?C + f t «) 

*4* = 

*** = S f' - < ) 

(b) Aerodynamics 

= -$ | [ h{Wt+e**Ha*s+*w}+*i(tA-$4,j} 

k >z = I [ r 3 { ru ♦« ft « 

+ ^ f e * :os ®>n + 6c s '" Sac)] 

k ,k *-% [ r i{ *s('-?)+itt%} + !,{*('- W‘)*s 
+ p, +e P T f 0 -f * p r 3 ) *< + 'P T r 0 3,- *v + «*&} 

+ r 2 {e^ 5 <- 3 t 2 *> 7 } - « c ] 
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S, = _iC, C« 4 2d 2 4 K + <*, 0 Ot + C 4 <^0 + W 
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t*i = ( c o +c t q,i)<x 
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Figure 1.- Schematic of the N-strips blade model considered; the hinge 
sequence is flap, lag, torsions; flap motion 8 positive up, 
lead-lag motion £ positive forward, torsion 6^ positive 
nose up. 
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Figure 4.- Root loci for a single strip (conventional) blade 
for o )0 = 5 and Xa 6 = 0 
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6.- Root loci for a single strip blade for wq = 0 and X Ag = .12c 
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Figure 7.- Root loci for a constant lift rotor (5 strips) blade for 
we = 0 X^ e = *12c 
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(c.g. coincidental with aerodynamic center) 



Figure 9.- Flap-lag roots for <u 0 = 0 and X Ag = .12c 
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Re X 

(b) Pitch-flap coupling Rg = -.3. 


Figure 12.- Flap-lag roots for u>e ■ 0 and * .12c. 
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(b) Pitch-lag coupling -.3. 

Figure 13.- Flap-lag roots for uq = 0 . nd = .12c. 
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Figure 14.- Flap-lag roots for u)q * 0 and X A = ,12c 
with tip Mach number Mp - .6 
(compressibility effect). 
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Figure 16.- Flap-lag roots for uig = 0, X a = ,12c, 
X c g = .06c, Rg = .-3 and structural 
lag damping g^ = .5. 
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